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Abstract 

The decay of dark matter particles which are coupled predominantly to charged 
leptons has been proposed as a possible origin of excess high-energy positrons 
and electrons observed by cosmic-ray telescopes PAMELA and Fermi LAT. Even 
though the dark matter itself is electrically neutral, the tree-level decay of dark 
matter into charged lepton pairs will generically induce radiative two-body decays 
of dark matter at the quantum level. Using an effective theory of leptophilic dark 
matter decay, we calculate the rates of radiative two-body decays for scalar and 
fermionic dark matter particles. Due to the absence of astrophysical sources of 
monochromatic gamma rays, the observation of a line in the diffuse gamma-ray 
spectrum would constitute a strong indication of a particle physics origin of these 
photons. We estimate the intensity of the gamma-ray line that may be present 
in the energy range of a few TeV if the dark matter decay interpretation of the 
leptonic cosmic-ray anomalies is correct and comment on observational prospects 
of present and future Imaging Cherenkov Telescopes, in particular the CTA. 



1 Introduction 



The existence of dark matter is now established beyond reasonable doubt by a variety 
of independent observations [TJ. These require the presence of substantial amounts of 
non-baryonic dark matter at vastly different scales ranging from individual galaxies to 
superclusters and filaments. Despite the overwhelming amount of gravitational evidence, 
however, no unambiguous evidence for non-gravitational dark matter interactions has 
been discovered to this day. Since a determination of the particle nature of the dark mat- 
ter from its gravitational interactions alone is impossible, searches for non-gravitational 
signatures of dark matter are of paramount importance. 

One of the principal approaches to its identification is the indirect detection of dark 
matter via searches for exotic components in the cosmic radiation produced by dark mat- 
ter interactions with Standard Model particles. For weakly interacting massive particles 
(WIMPs), it was pointed out some decades ago that the dark matter self-annihilation 
processes that can yield the correct thermal relic abundance might still occur today at 
a rate high enough to give rise to a flux of cosmic rays and photons that may be de- 
tectable by suitable telescopes. However, the self-annihilation of dark matter is not the 
only possible scenario for indirect dark matter detection. Namely, the dark matter might 
be unstable and decay into Standard Model particles, even though it must have a very 
long lifetime in order to survive from its production in the early Universe to the present 
day. If the dark matter is not perfectly stabilized by some unbroken symmetry, however, 
the possiblity exists that its decay products may leave visible traces in cosmic-ray fluxes. 
Indeed, in some well-motivated models the dark matter particles are not perfectly stable, 
but decay with cosmological lifetimes (see, e.g., ElElIllElinilZlIHlISl HD1 HH H21 US])- 

Among the possible indirect detection channels, measurements of cosmic-ray anti- 
matter are particularly sensitive to exotic contributions from dark matter. Interestingly, 
the PAMELA telescope has recently confirmed the existence of a dramatic rise in the 
positron fraction extending up to energies of at least 100 GeV [H], in stark contrast with 
expectations from conventional models of cosmic-ray production and propagation [13] . 
Furthermore, the Fermi Gamma-ray Space Telescope has observed a total flux of elec- 
trons and positrons that is harder than expected [HI [17], also indicating the possible 
presence of an additional component of charged cosmic-ray leptons [T8] . 

Various astrophysical explanations for these unexpected behaviors have been pro- 
posed [HI [20[ [21]. Arguably more exciting, however, is the possibility that the excess 
of positrons and electrons is due to the annihilation or decay of dark matter particles. 
However, measurements of cosmic-ray antiprotons, in particular measurements of the 
antiproton-to-proton ratio by PAMELA [22[ [23] , yield stringent constraints on the frac- 
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tion of dark matter decays or annihilations into hadronic final states. This has lead some 
authors to consider 'leptophilic' models of dark matter 0|211I251I2SII2ZII2SII2SII3DII3III32], 
where the dark matter is coupled predominantly or exclusively to charged leptons. In the 
following, we consider the possibility that the dark matter particles are indeed unstable, 
but decay leptonically with extremely long lifetimes. More precisely, the interpretation 
of the leptonic cosmic-ray anomalies observed by PAMELA and Fermi in terms of dark 
matter decay suggests a lifetime of the dark matter on the order of 10 26 seconds (see, e.g., 
[53]). This lifetime exceeds the age of the Universe by nine orders of magnitude, thus 
leaving the dark matter sufficiently stable on cosmological timescales. Nevertheless, due 
to the large amounts of dark matter in the Universe, even for such enormous lifetimes the 
resulting fluxes can be in the observable range. Indeed, strong constraints on decaying 
dark matter have been derived recently from gamma-ray observations of galaxy clusters 
and nearby galaxies [M] . 

In this paper, we examine some of the effects of leptophilic models of decaying dark 
matter which arise at next-to-leading order in perturbation theory and show that they 
can have relevance to indirect dark matter searches. We will not speculate on the pre- 
cise nature of the particle physics that could give rise to leptophilic dark matter decay. 
Instead, our approach will be to examine simple models where we assume effective inter- 
actions that describe the desired leptophilic coupling of dark matter particles to charged 
leptons. The salient point for us here is that even if one assumes an exclusive coupling 
of the dark matter to charged leptons at tree level, this behavior is only valid at leading 
order, while at next-to-leading order other particles, including photons and weak gauge 
bosons, will be produced. Indeed, these higher-order corrections have been analyzed in 
the past for the case of annihilating dark matter [351 EHJ [371 EH [391 SO] • It is well known 
that the higher-order corrections in the form of internal bremsstrahlung or from final- 
state radiation of weak gauge bosons can even dominate under certain conditions [361139] . 

The decay modes induced by higher-order corrections are usually suppressed by pow- 
ers of the couplings and possibly loop factors, as opposed to the leading-order decay 
modes. This means that the resulting decay products will be difficult to detect unless 
they possess some distinct features. Weak gauge bosons can be produced, for instance, via 
final-state radiation off the charged leptons [H]. By their subsequent hadronization, the 
massive gauge bosons will then generate hadronic particles, including antiprotons [40J. 
Therefore, every leptophilic dark matter model that aims to explain the leptonic cosmic- 
ray anomalies also serves as a source of antiprotons. 

In this work, however, we will focus on complementary constraints arising from a 
different decay channel induced by higher-order effects. Namely, we will study radiative 
two-body decays involving photons. These are particularly interesting, since they give rise 
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Figure 1: Tree- level diagrams contributing to the three-body decay ipr>M — ► i + i~N of 
fermionic dark matter, mediated by a heavy charged scalar S. Instead of the intermediate 
scalar £, the decay can also be mediated by a vector V. 

to monochromatic lines in the diffuse gamma-ray spectrum or in extragalactic sources. 
Such lines are of utmost importance because astrophysical processes generally gener- 
ate continuous gamma-ray spectra. Thus, the observation of a gamma-ray line would 
be a compelling signature of an underlying particle physics process. In some cases, a 
gamma-ray line can be produced already in tree level decays [12]. In the present work, 
we demonstrate that for leptophilic models of dark matter, the ratio between leading- 
order and next-to-leading-order decay modes can be large enough to produce a potentially 
observable gamma-ray line signal. 

This paper is organized as follows. In section [21 we discuss the production of 
monochromatic photons from radiative two-body decays induced at the one-loop level 
for fermionic dark matter particles in a simple leptophilic toy model. In section [31 we 
examine the corresponding case for a scalar dark matter particle. Next, in section [H we 
discuss observational constraints on gamma-ray lines in the GeV to TeV region and com- 
pare existing bounds with the expected signal from dark matter decay. We also comment 
on future observational prospects, in particular for the proposed Cherenkov Telescope 
Array. Finally, we present our conclusions in section [5j 

2 Radiative decay of fermionic dark matter 

We first regard the case that the particles comprising the dark matter are fermions 
which we denote by ipuM- We require that the dark matter decays with a large branching 
fraction into pairs of charged leptons in order to explain the excess of such leptons in 
high-energy cosmic rays, and we assume that this is the only channel in which the dark 
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Figure 2: Diagrams contributing at one loop to the radiative two-body decay ^dm — > 
jN, induced by a charged scalar £ (top row) and a vector particle V (bottom row), 
respectively. There are two additional diagrams in each case which differ only by the 
direction of the charge flow. 
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matter decays at leading order. If the dark matter carries spin 1/2 Lorentz invariance 
requires the decay to be (at least) a three-body decay involving a third, electrically 
neutral fermion N for angular momentum conservation. Thus, the decay i/jdm — > £ + £~N 
is the simplest one allowed by gauge and Lorentz invariance. Here, N could be a neutrino, 
a neutralino or a gravitino, for instance. The decays may be mediated by a virtual charged 
scalar particle S or by a charged vector boson V, with masses and my, respectively, 
which are assumed to be larger than the mass of the dark matter particle. We regard 
the two cases separately. 

In the case of an intermediate scalar, the effective Lagrangian that we use is given as 
the sum of a term coupling the dark matter, which we take to be a metastable Majorana 
fermion, to a charged lepton and a E particle, as well as a term coupling the neutral 
fermion to the E and a lepton field. We decompose the couplings into left- and right- 
handed components to allow for chiral couplings. Then the Lagrangian has the form 

Cf s = -^dm [\%Pl + \%Pr] IY^-N [\ L tN P L + \? N P R ] £ Et + h.c. , (2.1) 

where Pl = (1 — 7 5 )/2 and Pr = (1 + 7 5 )/2 are the left- and right-handed chirality 
projectors, respectively. The A-couplings can in general be complex. To obtain the re- 
quired cosmological lifetime for the dark matter, the couplings have to be super-weak 
or the mass of the mediator has to be super-heavy. The operators of the effective 
Lagrangian induce three-body decays of the dark matter into a pair of charged leptons 
and a neutral fermion at tree level, Vdm — > £ + £~N . The corresponding diagrams are 
shown in fig. [TJ 

In the case of a vector interaction, on the other hand, we assume an effective La- 
grangian of the form 

>Cff = -^dm7 m [X^Pl + A&PjJ I Vl - NY [X^ n Pl + X&Pr] £ V* + h.c. . (2.2) 

The case of mediation by a vector boson is more involved than the previous case of 
mediation by a scalar. In choosing a Lagrangian of this form, we assume that the essence 
of the gauge interaction giving rise to this Lagrangian is captured by the effective charged- 
vector interaction. In general, one expects neutral currents in association with the charged 
currents, which introduces a high degree of model dependence. For simplicity, we assume 
here that the decay is dominated by the charged-current interaction. 

1 For the case of spin 3/2, see refs. (43J [44] on radiative gravitino decay. 



6 



2.1 Decay widths 



In the following we examine decay modes of the dark matter at the tree- and one-loop 
level and summarize the relevant decay widthsj§ 

2.1.1 Tree-level decay: Vdm £ + £~N 

The leading-order decay induced by the effective Lagrangian is the three-body decay 
V'dm — > £ + £~N. If the dark matter decays in this way, it constitutes a possible explanation 
for the observed cosmic-ray anomalies under certain conditions [33]. We present the 
relevant expressions for the decay widths in the following. 

Mediation by a scalar. In the plausible limit rri£ <C rn^ DM <C mj, the partial decay 
width for the decay ^dm — > £ + £~N is given by (see app. Eland ref. [47] ) 

r(^ DM tl'N) = ^y 3 ^ {CfF^ml/m^J + Cf F 2 {m 2 N / m 2 ^)} . (2.3) 

The constants Cf, Cf are determined by the couplings as 

Cf = (|^| 2 + |A^| 2 ) (l^«v| 2 + l-^«v| 2 ) — V^ e (A^A^A^A^) 5 (2-4) 
Cf ee 2 v Re [(A^A^) 2 + (A* A^) 2 ] . (2.5) 

Here, r\ = ^ dm j]jv = ±1 depending on the CP eigenvalues of ipr>M and N. The kinematical 
functions, on the other hand, are given by 

Fx(x) = (1 - x 2 )(l + x 2 - 8x) - 12x 2 ln(x), (2.6) 
F 2 (x) = y/x[(l - x)(l + lOx + x 2 ) + 6x(l + x) ln(x)]. (2.7) 

In the hierarchical limit m^/m^ DM — > 0, the kinematical functions satisfy 

Fi(x) ~ 1, F 2 (x) ~ Jx for x -»■ 0, (2.8) 

whereas in the degenerate limit rriN/m^ DM — >■ 1, one gets 

~ -(1 - x) 5 , F 2 (x) ~ — (1 - x) 5 for x -> 1. (2.9) 
5 10 

In the limit m N <C m^, DM the decay rate (12. 3p corresponds to a lifetime 

— -*^(s)ra 5 (Tra) 4 - <-°> 

In the case where m^ DM and are quasi-degenerate, the decay rate scales approximately 
like (m^ DM - m N f. 



2 We have cross-checked the matrix elements for the three-body decays and the decay rates in the 
following sections by comparing them to the results from FcynArts |45| and FormCalc |46j . 
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Mediation by a vector. For the vector-mediated decay we find for the three-body 
decay rate in the limit mi <C rn^ DM <C m v 

1 AfnPi 

T(^ M -> trN) = 64(2?r)3 {CYF^Jm^J + C^F 2 (m%/ml DM )} , (2.11) 
where the functions F\ and F 2 are the same as in eqs. (12.6p . ( 12. 7p and 

Ci = (|A^,| 2 + l^l 2 ) (l-^«vl 2 + l-^«vl 2 ) + 2?? Re (A^A^-A^A^) , (2.12) 

Cl ee 2 V Re [(\^\%T + (A^) 2 ] = ( 2 - 13 ) 

For myv <C m D M anc ^ an analogous choice of parameters, the lifetime is smaller by a 
factor of four compared to eq. (12.10p . 



"x ^.(<£)(IM)\_~_)\ ,2,4) 
2.1.2 One- loop decay: Vdm - >• 

By combining the external charged lepton lines from the tree-level diagrams into a loop, 
we obtain diagrams contributing to the two-body decay -0dm — > 'jN (see fig. [2]). This 
decay mode will be suppressed with respect to the tree-level three-body decay by a 
loop factor and an additional power of the electromagnetic coupling. This is partially 
compensated by phase-space factors, however. More importantly, the two-body decay 
gives rise to monochromatic photons at an energy 

which can result in a distinct observational signature at gamma-ray telescopes, as will 
be discussed in some detail in section HI Interestingly, the experimental constraints on 
the parameters of decaying dark matter stemming from the non-observation of energetic 
gamma-ray lines could, despite the loop-suppression, be more stringent than the ones 
stemming from measurements of cosmic-ray electrons and positrons. 

Based on gauge invariance, and irrespective of whether the decay is mediated by a 
scalar or a vector particle, the matrix element for the sum of all diagrams contributing 
to the radiative two-body decay can be written in the following form, introducing an 
effective coupling B9> 

M = l -^u(k x ){P R - VNV ^ u P L )a^k 2ll e:u(p) 

= -^^HhKPR - VNV4, DM P L )hfu(p) , (2.16) 
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where = 7^]/2 and % DM , ?7tv are the CP eigenvalues of i/)dm and N, respectively. 
The partial decay width for *0dm — >• 7^ can then be easily calculated to be 



r(^ DM -+ IN) = %^m, DM ( 1 - ) . (2.17) 



The effective coupling g^^/i encodes all the information about the interaction between 
dark matter and the decay products. We give explicit expressions for this coupling in the 
following. 

Mediation by a scalar. We first examine the case of mediation by a charged scalar 
particle S (top row of fig. [2]). Assuming that CP is conserved in the interactions of ^dm 
and N, i.e., when the A-couplings are assumed to be real, the explicit form of the effective 
coupling gfj^ can be expressed as follows, 

+ (A&rA^ - v A&Aj) [v^m^um ~ K) - T) N m N K] } , (2.18) 

where the loop integrals /, I 2 and K are defined in app. lAlR The sum runs over all 
lepton flavors I G {e,/i, r}, for which Qi — Ci — 1. If multiple mediator particles (like 
left- and right-handed sleptons, £ = £l,£r) are present, one also has to sum over their 
contributions. In principle, there can also be contributions from quarks in the loop. 
Then, the sum runs over quarks and leptons with electric charge Q q and color charge 
C q = 3. However, tree-level decays into quarks can potentially lead to an overproduction 
of antiprotons if the relative size of the effective coupling to quarks compared to the 
coupling to leptons is too large. The requirement of avoiding antiproton overproduction 
then leads to the assumption of a leptophilic structure. For this reason, we assume 
throughout this work that the dark matter decays only into leptons at tree level. 

If the mass of the intermediate particle £ is much larger than the other masses, the 
loop integrals take on a very simple form. The effective coupling is then given approxi- 
mately by 

«5* - ^>i» f 1 " £^?{W« A » " ■ (219) 

For the concrete case where there is only one mediator S, which couples exclusively 



3 Note that the superscript '2' in the integral I 2 is an index, not a square. 
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to leptons, the decay rate reads 

r(^DM jN) 



m 



♦dm 



m 



N 



87r(647r 2 )^ m| 



m 



♦dm 



r]m N 

m V"DM 



X 



(A^yA^ 77 A^A^,) 



(2.20) 



For m^r mj, DM , this decay width corresponds to a partial lifetime 



tVdm-^JV - 7 x 1q29 s 



0.1 



[(A^-iyA^Aj)]" Wdm 



1 TeW m s \4 
10 15 GeW 



(2.2i; 



Mediation by a vector. In the case of mediation by a charged vector boson (bottom 
row of fig. [2]) we obtain the following expression for the effective coupling jH] , 

>M {(^DM^w^Vftf ~~ A^vA^,) [ij^ M m fcM (/ 2 - J - K) 



+ 77 JV m JV («/ - K)] + 2m,(^ DM A^Af^ - r] N \f N \^)J j , (2.22) 

where we encounter an additional loop integral J, which is defined in app. [A] Again, it 
is possible to include quarks by the replacement J2e ~~ J2f QfCf- However, as discussed 
above, this would not correspond to a leptophilic model. In the limit nig — > 0, m^ DM <C 
my, the above expression simplifies to 



v ^ 3e7] m^ pM / _ V m N 
9n t* 32^ m v \ m^ M 



(A^yA^, — rj A^A^,) 



(2.23) 



Thus, in the limit m^ <C m^ and m,/, DM <C my we obtain for the decay width 



9e 2 



♦dm 



m 



JV 



8tt (32tt 2 ) 



m? 



♦dm 



r/rriN 



m 



♦dm 



(2.24) 



For m,N m^, DM , this yields a partial lifetime 

0.1 



2 x 10 28 s 



1 TeV\ 5 / m v 
&MW r AM)]^^J vl^GeV 



(2.25) 



2.2 Intermediate scalar: intensity of the gamma-ray line 

The detectability of a loop-induced gamma-ray line will depend crucially on the ratio 
between the three-body decays at tree level and the two-body decays at the loop level. We 
examine the general expressions first and then evaluate them for some specific examples. 
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2.2.1 General expressions 



In the intermediate scalar case, the ratio between two- and three-body decay widths 
reads, neglecting the charged lepton masses, 

r(^DM^7^) _ 3a cm [E< ( - y xf N x%)] 2 (i - aQ 3 (i - vV^) 2 



where x = m^/m^ DM and the kinematical functions F\ and F 2 were defined in eqs. (12.61) . 
(12. 7p . This general expression can be used to study the intensity of the one-loop in- 
duced gamma-ray line in different scenarios. The numerical value of the prefactor is 
3a cm /(87r) ~ 1/1148. 

In general, the fraction depends on the chiral and flavor structure of the couplings, 
the mass ratio Ttjv / wi^ DM of the decay product and the dark matter particle, and the 
relative CP parities r\ = ±1 of N and ^dm- For many practical purposes, it turns out 
that the dependence on the couplings \m/if> an d on kinematics, i.e. on x = m^/m^ DM , 
can be factored according to 



57T 



x R (X m , X ei p, X eN , \m) x S^mpf/m^^g) . (2.27) 



In this parametrization R% captures the model-dependence, whereas S v is determined 
entirely by kinematics. 

It is interesting to consider the two limiting cases of hierachical masses, rriN /m^ UM — > 
0, and degenerate masses, /m^ DM — > 1. In the hierarchical limit, and assuming for 
simplicity real couplings, one explicitly obtains S^ lcr = 1 and 

2 



J2e [\^e$ + A f$) { A w + A w) _ V A e4> A fN A i% A w] 



p£,hier/\Z \L \R \R \ L^f •l"lN' y lj>)\ /n 98 \ 

U v {A m , A m , Aty) — , \R2\(\L2 , \R2\ „ \L \L \R \R 1 " 



For generic couplings, i?^' hlcr is roughly of order one, unless for some special cases where 
cancellations or chirality suppressions occur. It follows that in the hierarchical limit the 
two-body decays into 7iV are typically suppressed roughly by a factor 10~ 3 compared 
to the tree-level decays into £ + £~N. In the next subsection we will examine the model- 
dependent factor R^j for some specific cases. 

On the other hand, in the degenerate limit /fn^ DM —> 1, and again assuming real 
couplings, one finds 

-R^' deg (A^v,A^,, \f N , Xf^) = (2.29) 

[J2i { A £N A £iP ~ V A £N A fi>)] 
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which is also roughly of order one for generic couplings, i?J- deg ~ Oil), and 

^(-/™,»)-{ 20/(1 _ 5 4 2 KdJ2 iv.i. ™ 

Thus, for the case rj = +1, i.e. when i/jbm and N have the same CP parities, we again 
find a typical suppression factor of the order of 10~ 3 for the two-body relative to the 
tree-level decay rate, as in the hierarchical case. Interestingly, however, when ipBU and 
TV have opposite CP parities, r\ = — 1, the two-body rate can be enhanced significantly 
even for a relatively mild degeneracy, as is shown in fig. [31 This enhancement is due to 
the fact that the decay rate r(^ DM —> jN) is proportional to [m N — m^ DM ) 3 , whereas 
the decay into leptons is suppressed like (mjv — m^ DM ) 5 [48 J. Most interestingly, due to 
this enhancement the decay rate into '-/N can be rather large in some cases, yielding 
potentially very intense gamma-ray lines. 

As a side remark, we note that in addition to the decay channel ^dm — > 7^ into 
photons, there can exist a decay mode ipr>M — > Z°N into Z-bosons, which can naively 
be expected to be of similar size. Thus, for situations where the gamma-ray line signal 
is strongly enhanced, an equally enhanced decay into Z-bosons can yield additional 
constraints from the antiproton flux produced by the subsequent fragmentation of the 
Z-bosons. We leave a more detailed discussion for the future 1491. 



For concreteness, we consider the case of purely chiral, say left-handed, couplings, 
and that only one mediator species is present. Then the fraction of decay rates is given 
by 

{1 for rriN — > 0, rj — ±1 
V2 for^^ DM ,, = +l . 
„ 20 for TTi ^ Tfi 71- 1 

(2.31) 

In this case the model-dependent factor for the hierarchical and degenerate regimes 
coincides, _R s ' hier = ^ s > de s = _R chir ; and is furthermore independent of rj. Explicitly, one 
has 



The result for purely right-handed couplings is analogous. For a generic choice of cou- 
plings, one expects R chn ~ 0(1). Note that R chlv < Ng, where Ng is the number of flavors 
participating in the decay. 

For example, consider two particular cases for the flavor composition of the lepton 
pairs produced in the decay: 
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Figure 3: Ratio of the decay rates r(^DM — > jN)/ Xlf T(^dm — > i + i~N) when the 
decay is mediated by a scalar. The four cases correspond to single-flavor decay (red) and 
democratic decay into all flavors (blue), as well as iPbm/N having the same CP parity 
(dashed, r] = +1 ) or opposite CP parity (solid, r\ — —1 ). See eq. (12 .311) . 

(A) Decay into a single lepton flavor: 

(B) Flavor-democratic decay into e + e~, t + t~. 

Then one has R chlT = 1 in case (A) and R chn = 3 in case (B). 

The dependence of the ratio of decay rates on the mass ratio /m^, DM is shown in 
fig. [3] for the two cases (A) and (B), and for rj = ±1. This dependence is in fact a rather 
generic feature, which is independent of the details of the couplings. We emphasize again 
that, in the case when V'dm and have opposite CP parities, t] = — 1, even a rather 
mild degeneracy between and m^ DM can lead to a considerable enhancement of the 
gamma-ray line signal relative to the electron/positron flux. 

2.2.2 Examples 

Now, we will discuss the branching ratio into gamma-ray lines for several specific sce- 
narios. Namely, we consider the case where N corresponds to left-handed neutrinos ul, 
as well as the scenario of kinetically mixed hidden U(l) gauginos, where N corresponds 
to a neutralino. 
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Decay into left-handed neutrinos. 

As a basic example for the scalar-mediated decay described in the previous subsec- 
tion, we consider the case where the neutral fermion is a left-handed neutrino, N = vl- 
Then one can set \\ v = and m v = 0. 

From eq. ( I2.26|) it directly follows, that in the limit ra t <C ?7i^ DM <C the ratio 
reads 

r(^DM ->• 7v) 3a em [Y^e A iv A hiA 



,2 



(2.33) 



As long as only one virtual scalar particle is relevant for the three-body decay, the last 
factor in this expression is bounded from above by the number of lepton flavors Ni 
that contribute to the decay. For N(, = 3, this corresponds to a branching ratio into 
monochromatic photons smaller than 3 x 10 -3 . 

Hidden-gaugino dark matter. 

In supersymmetric scenarios with an extra U(l)x gauge group in the hidden sector, 
which kinetically mixes with the Standard Model U(1)y, the particles ^dm and N could 
be associated with the hidden gaugino and the lightest MSSM neutralino, respectively. If 
the kinetic mixing parameter 9 is extremely small, the hidden gaugino could constitute 
decaying dark matter [9]. For a bino-like lightest neutralino and E = £l being a left- 
handed slepton, the couplings are approximately given by 



* ^f £ L 9, (2.34) 



X L m ~ ±=Y t L , (2.35) 

where 9 ~ 10~ 24 is the mixing angle of hidden gaugino and bino, fixed by the requirement 
of a lifetime of the order of 10 26 s, and Y e L = +1. 

For a hidden gaugino that decays into a bino-like neutralino, the decay rates are 
given in app. A of ref. [9]. One can also obtain these rates using the expressions derived 
above. In particular, we have to sum over two 'mediators' £ for each flavor: 

m v — Tt- \ L — 9' v l q \ L — 9' v L \ R ~ n \ R ~ n 

• ^ _ l L- A ty — -jl 1 1 v i A IN — ~j2*l ' A lip — u ' A IN — u 

. E = £ R : AJ = -faYfy xf N = -far*, \% ~ 0, \ L m ~ . 

In addition, there are corresponding contributions from (s) quarks. We assume mj ^> 
m^ DM , and neglect the mixing of for simplicity. Note that there are additional 
contributions from chargino loops [IB] that are suppressed by the fourth power of the 
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inverse chargino mass. We assume that the squarks are much heavier than the sleptons, 
and we furthermore assume that all slepton masses are degenerate. Finally, if we take 
the limit mt — > 0, we get 



3 - -2 



87r(327r 2 r 16mi \^ m^J V ™</> D 
x [3(1 -4r/)] 2 . (2.36) 

If the bino and the hidden gaugino have the same CP eigenvalue, one has r\ = +1, 
otherwise 77 = —1. 

The three-body decay rate can be obtained from eq. (12.31) . which can be easily gen- 
eralized to also account for neutrinos and quarks in the final state. Note that in general 
one has to add the matrix elements for the decays mediated by £ = II and X = and 
compute the decay rate from the square of the summed matrix elements. However, it 
turns out that all 'interference '-terms are suppressed by the bino-higgsino mixing^] which 
we neglect here. Thus, it is possible to add the decay rates directly. Note that there is 
an additional contribution from a Z° on the intermediate line [17], which is subdomi- 
nant for the parameter range considered in ref. [9]. Therefore, we also neglect it here for 
simplicity. We assume, as above, degenerate sleptons. The decay rate summed over three 
generations of charged leptons and neutrinos is thus (assuming m-£> ~ mj) 

£r(^ DM -+ ON) - x 3 x 18 x (F, + 2 V F 2 ) . (2.37) 



In the hierarchical limit <C m^ DM , the kinematical factor approaches unity, F\ + 
2r\F 2 — > 1. In the degenerate limit m N — > m^ OM) one finds Fi + 2r]F 2 —> (2 + 77) (1 — 



The ratio of decays into 7iV to the decays into charged leptons is thus given by 

r(^ DM -+ jN) 3a cm [3(1 -4?7)] 2 (1 - xf (1 - Wi) 2 

Et^(thM^£ + £-N) 8vr 51 F l (x) + 2r ) F 2 {x) ' { ' ' 

* „ -v ' 

= Rr) —Sr] 

with x = m1/m^ DM . For the hidden gaugino, the decays -^dm — > Z°N and ^dm — >• h°N 
can also be important since they occur at tree level. Their rates are given in eqs. (A.l) 
and (A.3) of ref . 0. 

4 The reason is the following: If we consider a pure bino-slepton-lepton interaction, the slepton II 
couples only to left-handed leptons and the slepton In only to right-handed ones. Thus both channels 
are 'orthogonal' in the limit where neutralino and slepton mixing are neglected. 
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According to eq. (12. 38 p . the model-dependent factor R v is here given by R + = 1.6 for 
r) = +1 and i?_ = 4.4 for 77 = —1, respectively, and hence of order one. The kinematical 
factor S v is precisely of the form that was discussed in section 12.2.11 where we found 
that the two-body decay rate may gain significantly in importance relative to the three- 
body decay rate if the masses of the hidden gaugino ipT>u and the neutralino N are 
near-degenerate and the two particles have opposite CP parities. 

2.3 Intermediate vector: intensity of the gamma-ray line 

In the case of mediation by a vector, the ratio between two- and three-body decay rates 
is 

T^dm Ny) _ 27a cm [£, {Xf N X^ - r, \f N \%)] 2 (I - xf{l - r^/xf 



Ee T ^DM -+ £ + i-N) ~ 8tt T.t(PXF 1 {x) + G\ F 2 (x)) 

where x = m 2 N / m^^. As before, it is useful to consider the hierarchical limit 
m7v/m^, DM — > 0, and the degenerate limit rriN/m^ DM — > 1, for which it is possible to 
capture the dependence on the couplings in a factor BY and on kinematics in a model- 
independent factor S v , 

x R n ( A W' X U %• x ^("W™Vdm) • (2-40) 



7T 



The kinematical factors S v are identical to the case of scalar mediation, see eq. ( 12.301) . 
For the model-dependent factors, one finds 



( A w A ej, V A fN A ?j>)] 

n v 1%' A £il>i A Mi A £^) - ^ r/vL2 , \R2\(\L2 , \R2\ i 9t) \L \L \R \R 1 ' ^■^ L > 

1^1 \\ A tip + A eip)\ A w + a in) + lr l A e4< A iN A i^ A iN\ 

in the hierachical case, and 



R^' des (\f N , \f N , — (2.42) 

[J2i { A £N A fip ~ V A fN A hjj)] 



EI\L2\L2 i \R2\R2\ , 2_(\L2\R2 _i_ \R2\L2\ _i_ \_R \L \_R 

« \ A tip A W ~T A li> A lN) 2+r)\ A li> A tN A lil) A lN) 2+r, A (.ip A Hip A tN A iN 

in the degenerate case. For purely chiral, e.g. left-handed, couplings one finds that 
i?^' hlcr = i?^ ,dcg = i? chir coincides with the expression ( 12.32ft for the scalar case, as 
does the kinematical factor. For a generic set of couplings, Rl is roughly of order one. 

Note that the prefactor of the ratio of decay rates, eq. (12.401) . for mediation by a 
vector is larger by a factor of nine compared to mediation by a scalar, eq. (12.27]) . Thus, 
in the hierarchical case m^-/m^ DM — > as well as in the degenerate case with t] = +1 
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Figure 4: Ratio of the decay rates r(-0DM — > J2iT(ipDM — > £ + P N) for decay 

mediated by a heavy vector. Otherwise, the four cases are identical to the ones shown in 
fig. 13 

one finds a ratio between two-body and tree-level decay of the order of 10 -2 , one order of 
magnitude larger than for the scalar case. In addition, when rj — — 1 the gamma-ray line 
is further enhanced for m N /m^ DM — > 1 by the kinematic effect discussed in section [2.2. II 
The ratio for some specific examples is shown in fig. HJ 

Therefore, there are scenarios with dark matter decay mediated by heavy vectors 
where a gamma-ray line can be fairly intense, despite being loop-suppressed, while at 
the same time being in agreement with the electron/positron measurements. 

3 Radiative decay of scalar dark matter 

We now consider the case that the dark matter particle is a (pseudo-) scalar which we 
denote by 0dm- In this case, the symmetries allow for the decay into a pair of charged 
leptons at tree level, 0dm _ > We describe this by an effective Lagrangian that 

describes a direct interaction between dark matter and charged leptons, 

£cff = -I [X^Pl + \%Pr] i 0dm + h.c. . (3.1) 

If the dark matter particle is a parity eigenstate, one has = Xf^ for a scalar and 
A£ = — Xf, for a pseudo-scalar. 
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Figure 5: Tree- level decay of scalar dark matter. 




Figure 6: Diagrams contributing to the two-body decay of scalar dark matter into two 
photons at the one-loop level. 
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3.1 The decay </> DM £ + £~ 

The effective Lagrangian (13.11) will give rise to the tree-level decay shown in fig. [5j The 
corresponding decay width is 



r (0dm 



\Mf 



Am 2 , 



2 ' 



(3.2) 



where m^ DM and mg are the mass of the dark matter and the charged leptons, respectively, 
and the amplitude is given by 

\M\ 2 = ml DM (\\$/+\\fX)-2m}\\^ + \?/ ■ (3.3) 

Thus, in the case of equal left- and right-handed couplings, X\^ = Xfi = A^, one gets 

3/2 

«-) = ^oJl-^) ■ (3.4) 



r (0dm ->■ 



8tt 



??? 



</>DM 



For mi <C tti^ dm this corresponds to a lifetime 



10 -26X2 / 1TeV \ 

^ s2xl0 " B i-) 



(3.5) 



3.2 The decay 0dm — >• 77 



By combining the external lepton lines into a loop, decays into two monochromatic 
photons radiated off the charged lepton loop are induced at the quantum level (see 
fig. ED- 

For equal left- and right-handed couplings, Xf^ = Xf^ = A^, and in the limit me <C 
m 4>DM ( see app. [B]and [501 I5T] ). 



-»■ 77) 



0DM 



.2 \ 2 



167T \ 167T 2 

where for the relevant limit r ^> 1 one has 

A/(r) ~ 



I |2-i(ln(4r)-t7r; 



(3.6) 



(3.7) 



Thus, when taking only one lepton species into account, we obtain for the ratio between 
the decay into photons and charged leptons 



'DM 



->■ 77) , "em "1 



'DM 



27r 2 m\ 

0DM 



2 - i (ln(4r,) - ztt) 2 



10 



-9 



m e 



2 /l TeVV 



106MeV; Wdm/ 
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(3.8) 



We see that for scalar dark matter, the decay into two photons is highly suppressed by 
the factor ■mj/m'^^ compared to the decay into a pair of charged leptons. In addition 
to this helicity-suppression factor there appears a factor a^ m /7i 2 as opposed to a cm /n 
for fermionic dark matter, since the loop contains two photon vertices, and both the 
tree-level and one-loop decays are two-body decays. The same suppression factors occur 
for pseudo-scalar dark matter and for the decay into massive gauge bosons. Thus, there 
appears to be no hope of detecting a gamma-ray line in this case. For more general 
expressions for the decay rates, see app. [Bj 



4 Observational constraints 



The observation of a cosmic gamma-ray line at TeV energies would be a strong hint 
for the dark matter interpretation of the PAMELA/Fermi LAT e ± anomalies. On the 
other hand, the non-observation of gamma-ray lines can be used to constrain the above 
leptophilic models, which induce these lines at one loop, as discussed above. The gamma- 
ray lines that originate from dark matter decay inside the Milky Way halo could be 
observed in the isotropic diffuse gamma-ray flux. Furthermore, lines may be observable 
in the flux from nearby galaxies and galaxy clusters. 

At intermediate energies, satellite instruments such as Fermi LAT are a very sensitive 
probe for gamma-ray lines in the Galactic flux. At higher energies, Imaging Atmospheric 
Cherenkov Telescopes (IACTs) provide important information. For the future, the pro- 
posed Cherenkov Telescope Array (CTA) is expected to improve the flux sensitivity of 
current IACTs (MAGIC, H.E.S.S., VERITAS) by an order of magnitude. We put some 
emphasis on IACTs, since these instruments are capable of probing the high energy 
ranges relevant to the dark matter interpretation of PAMELA/Fermi LAT. 



4.1 Fermi LAT 

The flux of monochromatic gamma rays from the decay of dark matter in the Milky 
Way halo is given by a line-of-sight integral over the dark matter distribution [52] . This 
component of the gamma-ray flux is explicitly given by 



dJZ° _ T(^dm -yN) 



dE 4irm 



S(E,-E) J dtp^if), (4.i; 



4>DM 

l.o.s 



where T(^dm — > 7^0 denotes the partial decay width of dark matter particles for two- 
body decays involving a photon and a neutral particle N. When the neutral particle is 
massless, we will write v instead of iV in the following. Furthermore, fn^ DM is the mass 
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of the dark matter particle, E 1 is the energy of the produced gamma-ray line as given 
by eq. ( 12 . 1 5j) . while p^^f is the Milky Way's dark matter halo density profile. We adopt 
the NFW profile here, which has the form 

= . * . V2 , (4.2) 
r/r c (1 + r/r c ) 

with the parameters p c = 0.35 GeV cm -3 and r c = 20 kpc [5*31 151] . leading to a local dark 
matter density of 0.4 GeV/cm 3 [55]. The gamma-ray flux from dark matter decay inside 
the Galactic halo has only a mild angular dependence and can be considered as isotropic 
for our purposes (for details on anisotropies in the Galactic gamma-ray flux from dark 
matter decay, see refs. [521 EE])- The extragalactic contribution stemming from the decay 
of dark matter at cosmological distances is generally fainter than the Galactic flux, and 
we will neglect this component here. 

The Fermi LAT collaboration has conducted a negative search for Galactic gamma- 
ray lines in the diffuse flux in the energy range from 30 to 200 GeV [51] • For the halo 
profile (14. 2p . we plot the resulting 2a limits on the partial decay width corresponding to 
V'dm — > l v in fig. Most interestingly, the Fermi LAT observations can constrain the 
dark matter decay into photons at the one-loop level if the total dark matter lifetime 
is of the order 10 26 seconds. Thus, the Fermi LAT bounds on gamma-ray lines can be 
relevant for dark matter scenarios with m^ DM ~ 300 - 400 GeV, which can provide a 
possible explanation for the rise in the positron fraction observed by PAMELA (see, e.g. 
[57]). 



4.2 Imaging Atmospheric Cherenkov Telescopes 

IACTs are important tools to constrain scenarios with dark matter masses in the multi- 
TeV range. One property of these instruments is that the atmospheric showers induced 
by cosmic-ray electrons or gamma rays cannot be distinguished easily, since both particle 
species initiate similar electro-magnetic cascades in the atmosphere. The large cosmic- 
ray electron flux hence comprises an irreducible background for high energy gamma-ray 
observations. Since the electron background is expected to be very isotropic — in contrast 
to the gamma rays — it can be removed by calculating differences between fluxes that 
are observed in different neighboring regions of the sky. As a result, IACTs are best 
suited to observe localized sources, whereas diffuse signals such as those resulting from 
dark matter decay are more difficult to discern from the background unless they exhibit 
sharp spectral features. Constraints on the gamma-rays from decaying dark matter can 
be derived in two different ways. First, one can observe point-like sources like M31. 
Second, by using the observed electron+gamma-ray flux (potentially also contaminated 
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by unrejected protons), one can derive upper limits on the Galactic halo signal from 
dark matter decay. If the statistics are good enough, one could even hope to see spectral 
features in the electron+gamma-ray flux, or translate their non-observation into bounds 
on the corresponding dark matter decay width. This will be described in the context of 
the CTA below. 

The HEGRA collaboration has published constraints on the gamma-ray line flux 
from M31 [58]. These bounds can be converted into 99% C.L. limits on the decay width 
of dark matter into gamma-ray lines. HEGRA observed a region with an opening-angle 
of # bs = 0.105°, corresponding to the inner 1.4 kpc region of M31. The expected flux of 
gamma rays from dark matter decay from M31 can be derived as follows. We define 8 
to be the angle between the line-of-sight and the ray that passes through our position 
and the center of M31. Each angle 9 then corresponds to an 'impact parameter' R. If 
D is the distance to the target (D = 770 kpc in case of M31), we have R ~ DO. The 
gamma-ray flux from dark matter decay in M31 within the opening angle # bs is then 

M31 e ° hs °° 

% L = F( f ^ 7A ° 6(E y -E)2n [ dOsine [ ds f%g(V* + ff) , (4.3) 
dE 47rm DM J J 

-oo 

where the first integral is over the solid-angle, whereas the second integral is over the 
line-of-sight. For the dark matter density profile of M31 we adopt the NFW profile 
with values given in ref. [59], p c = 2.0 GeV/ cm 3 and r c = 8.31 kpc. The other profiles 
from ref. [59] lead to similar constraints. The signal from decaying dark matter has a 
relatively large angular extent due to the linear dependence on the halo profile, and can 
leak into the off-region which is used to estimate the background fluxes of the I ACT. The 
details of this effect depend on the details of the adopted off-region and are different for 
each observation. Here and below, we incorporate this effect simply by subtracting from 
eq. (14. 3 p a flux corresponding to the dark matter- induced flux emitted at 9 = 2 ^obs- This 
should lead to correct bounds within a factor of two. Our results are shown in fig. [71 

Upper limits on the gamma-ray flux from the Perseus galaxy cluster were presented 
by the MAGIC collaboration in ref. [60]. For the density profile of the Perseus cluster we 
take the NFW profile with r c = 384 kpc and p c = 0.04 GeV cm -3 , the obsevational angle 
is 6* bs = 0.15°, and the distance to the Perseus cluster is 78 Mpc. The resulting 95% 
C.L. bounds are shown in fig. [7]§| Since the energy threshold of the MAGIC telescope is 
very low, we can constrain gamma-ray lines with energies down to 100 GeV. 

The H.E.S.S. collaboration has published measurements of the electron flux at TeV 
energies [611 [62] . The measured electron flux may be contaminated with diffuse gamma 



3 We take the limits corresponding to T = —2.5 from tab. 4 of ref. [60] 
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— Fermi Halo -- HEGRA M31 
H.E.S.S. tr -■- MAGIC Perseus 



M DM [GeV] 

Figure 7: Lower bounds on the inverse decay width of dark matter decaying into gamma- 
ray lines via t/>dm — ► l v are shown as black lines. The bounds on this decay channel come 
from line searches in M31 by HEGRA, from line searches in the diffuse flux by Fermi 
LAT and from observations of the Perseus cluster by MAGIC. Further bounds can be 
derived from the (7+)e~ observations of H.E.S.S. Our estimates of the reach of the future 
CTA in measurements of the flux from M31 or spectral variations in the diffuse 7 + e~ 
flux are shown as red lines. 



# 


Channel 


V 


m DM [GeV] 


£ 7 [GeV] 


rriN [GeV] 


^ i+i-N N 


m DM r sTeV ] 


1 


ele+N 


-1 


1000 


170 


812.4 


2.5 x 10 26 


2.47 x 10 27 


2 


e- L e + L N 


+ 1 


500 


170 


282.8 


5 x 10 26 


6.51 x 10 29 


3 


eleiN 


-1 


400 


170 


154.9 


6.3 x 10 26 


4.83 x 10 28 


4 




-1 


100000 


5000 


94868 


4.5 x 10 24 


2.87 x 10 26 


5 




+1 


15000 


5000 


8660 


3 x 10 25 


1.18 x 10 30 


6 




-1 


15000 


5000 


8660 


3 x 10 25 


3.06 x 10 28 


7 




+1 


15000 


5000 


8660 


3 x 10 25 


3.42 x 10 34 



Table 1: Benchmark scenarios. In the first three cases, the three-body decay produces 
only electrons. In the last four cases, the three-body decay produces muons. The gamma- 
ray line intensity of these scenarios is illustrated in fig. O 
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ip DM -^l + £ v (intermediate scalar) 

10 30 r , , , — , — . .— 




M ml [GeV] 



Figure 8: Same as fig. [7J In addition, the orange and gray shaded regions show the 
parts of the parameter space that are relevant for the dark matter explanation of the 
PAMEL A/Fermi e anomalies with the flavor- democratic decay channel ^dm — > l + l~v. 
The intermediate particle is assumed to be a scalar, in which case the branching ratio 
into monochromatic photons can be as large as BR(^rjM ~> 7^) — 3 x 3a om / (8ir), which 
we assume here. 
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ip DM -^-£ + £ v (intermediate vector) 

10 30 r , , , — , — , ,— 




M DM [GeV] 



Figure 9: Same as fig. El but assuming that the intermediate particle is a vector, in which 
case the branching ratio into monochromatic photons can be as large as BR(^dm — * 
yv) ~ 3 x 27a cm /(87r), which we assume here. 




M DM [GeV] 



Figure 10: Same as fig. El but for decay into -0dm 1^ ~yTv. The intermediate particle 
is assumed to be a scalar, leading to a maximal branching ratio of BR(^dm — > 7^) — 
3a e m/(8vr), which we assume here. 



25 




Mm I GeV l 



Figure 11: Same as fig. El but for decay into -0dm — > /i + /i~zA The intermediate particle 
is assumed to be a vector, leading to a maximal branching ratio of BR^dm — > 7^) — 
27a em /(87r), which we assume here. 



rays by no more than « 50% [01]. This fact allows to translate the electron flux into upper 
bounds on gamma-ray lines from dark matter decay in the Galactic halo. For energies 
above 1 TeV, we derived 2<r-bounds from the fluxes shown in fig. 3 of ref . [01] • For energies 
below 1 TeV, where the H.E.S.S. results overlap with the Fermi LAT measurements of 
the electron flux, la-upper limits on the amount of diffuse gamma rays were derived 
by comparing the H.E.S.S. and the Fermi LAT electron fluxes in ref. [03]- These upper 
limits can also be used as bounds on gamma-ray lines. Our results are shown in fig. [7J 



Prospects for the CTA. We will now briefly discuss observational prospects for 
the future Cherenkov Telescope Array (CTA, see ref. [01] for a recent discussion). The 
expected 2<7-limit from M31 that the CTA could produce can be roughly estimated by 

/ tM31\ < max(2y/AU,3.1) 

Wdm /on rsj rp A 5 K^-^J 

J- ^eff 

where N OQ denotes the number of measured events in the on-region, T is the measurement 
time, A e g denotes the effective area of the instrument (we take A e g m 2 km 2 at 5 TeV, 
and let it scale with the energy as in ref. [OS], fig. 17a), and (J^^on is the gamma-ray 
flux from M31 averaged over the on-region. As on-region, we take a circle with 1.0° radius 
around the center of M31, and for the off- region we assume that the solid-angle of the 
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10 2 10 3 10 4 

£ 7 [GeV] 

Figure 12: Same as fig. [7J but with different scaling of the axes to allow for no n- vanishing 
m^. The black squares and red dots show the predictions for the different benchmark 
scenarios summarized in tab. [TJ Black squares correspond to scenarios with rj — +1, 
while red dots correspond to rj = — 1. The last benchmark point in tab. [T]lies outside of 
the shown parameter region. 
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off-region is much larger than the on- region, Q g ^> f2 on . If we assume that only the 
background is observed, and no signal is coming from M31, N on can be estimated by 

N on = n on TA cS (J e - + e r J p ) . (4.5) 

Here, J e - and J p denote the cosmic-ray electron and proton fluxes, respectively, and e r 
is the rejection factor of protons. Fluxes have to be integrated over an energy range that 
corresponds to the energy resolution of the detector (around 10%, taken from ref. [64J, 
fig. 23, scenario E). For the cosmic-ray electron and proton fluxes at high energies we 
take 3 

^ = 1.17 x 10- 11 ( ' cm" 2 s" 1 sr" 1 GeV" 1 , (4.6) 
dE \ TeV / 

from ref. [61], and 

^ = 8.73 x 1(T 9 ( —) 2 7 cm" 2 s^ 1 sr" 1 GeV" 1 , (4.7) 
dE \ TeV / 

from ref. [65], respectively, in agreement with the cosmic-ray measurements. Furthermore, 
at energies below 1 TeV, the electron flux becomes somewhat harder with a spectral index 
of ~ —3.1, and we replace eq. (14. 6 j) with a flux fitting the results of ref. [16] in this energy 
regime. Taking into account also other cosmic-ray species beside the protons would only 
have minor impact on our results. The proton rejection factor is set to e r ~ 10~ 2 [64"t [65]. 
and we use as observational time of M31 T = 20 h. Our estimates for the limits that the 
CTA could produce from M31 observations in the future are shown in fig. [7] by the lower 
dashed line. They are almost two orders of magnitude better than the limits derived 
from the HEGRA observation. This is mainly due to the increased effective area of CTA, 
but also due to the larger on-region that we adopted in our estimates. For decaying dark 
matter it is not optimal to search for point-source signals, as was done in the HEGRA 
analysis, for example. Using a larger on-region typically leads to better results. 

Instead of using the spatial variations of the observed cosmic-ray flux to derive con- 
straints on dark matter decay in extragalactic sources, one can also derive constraints 
from the non-observation of spectral line features in the diffuse flux, which could come 
from dark matter decaying into gamma-ray lines in the Galactic halo. In this case it 
is best to consider data from large fractions of the sky, to maximize the statistics. The 
expected "halo" -bound that the CTA will presumably reach then follows from 



(Jd -° )sky " YX^ ' (4 ' 8) 

where (J^°)aky denotes the gamma-ray flux coming from dark matter decaying in our 
Galactic halo, averaged over all angles. We assume that the data will be good enough to 
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estimate the background by fitting a power-law to the observed flux at energies close to 
the line, similar to the analysis in ref. [51], and we neglect the statistical uncertainties 
in the background estimate^ In eq. (14. 8\\ . N is the total number of observed events, 
including electrons, gamma-rays and protons that pass the cuts. The region Q is taken 
to be as large as possible to maximize the statistics (we assume Q = 7r(3°) 2 ), and as 
observational time we take T = 1000 h. As above, we integrate over energy bands which 
correspond to the anticipated energy resolution of CTA. Our resulting estimates for the 
bounds that CTA could obtain observing the diffuse flux are shown in fig. [7J by the 
upper dashed line. As can be seen from this figure, the bounds on gamma-ray lines from 
dark matter decay that can be put by looking at spectral variations in the observed 
diffuse fluxes can be even stronger than the ones that can be derived from flux limits on 
point-like sources like M31. 

4.3 Discussion 

The case m^v — > 0. We first discuss the case where dark matter decays into a photon 
and a massless particle. In fig. [7J we present a collection of the lower bounds on the 
inverse decay width for two-body decays into a monochromatic photon and a massless 
particle as determined by the methods described in the previous subsection. For dark 
matter masses between 100 and 400 GeV, the line searches in the diffuse Galactic flux 
by the Fermi LAT constitute the strongest constraints. At higher energies, Cherenkov 
telescopes provide important information. As far as constraints from particular sources 
are concerned, we show the constraints from HEGRA observations of M31 and MAGIC 
observations of the Perseus cluster. We also plot the constraints from the diffuse electron 
flux observed by H.E.S.S. Lastly, we show our estimates for the reach of the future CTA 
which could improve current limits by almost two orders of magnitude at energies above 
a few hundred GeV. 

In fig. [8] we show the same constraints together with shaded regions indicating the 
part of the parameter space relevant to PAMELA and Fermi for the gamma-ray lines 
induced by the decay ipBM l + l~v. The orange regions correspond to the fit to the 
positron fraction as measured by PAMELA, whereas the dark gray regions correspond 
to the fit to the total flux as measured by Fermi LAT. In both cases, the lighter 
shades indicate the 5a confidence level around the best-fit point, while the darker shades 
indicates the 3a confidence level. We only regard the data points above 10 GeV, which 
are not significantly affected by solar modulation. For the background fluxes of secondary 

6 Note that this is different from our treatment of the H.E.S.S. electron flux, where we only required 
that the predicted line signal is below the observed fluxes, without any attempt to subtract a power-law 
background. 
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electrons and positrons, we assume the 'model 0' backgrounds [18] as parametrized in [33] . 
In the energy range of interest, we assume that the primary electron flux is given by a 
simple power law. At each point in the (todm, TbM)-plane, we then allow the power- 
law index of the primary electron flux to vary between —3.0 and —3.3, whereas the 
normalization is fitted to the data. We find that the relevant parameter space is not 
constrained by current instruments, but could be constrained by CTA in the future. 

The same plot is shown in fig. [9J but assuming that the decay is mediated by 
an intermediate vector particle, in which case the branching ratio can be as large as 
3 x 27a cm /(87r), which we assume in the figure. This is about an order of magnitude 
larger than in the case of mediation by a scalar, and one can see that in this case the 
CTA can indeed constrain a significant part of the parameter space relevant to the dark 
matter interpretations of PAMELA and Fermi. Analogously, in figs. [101 [TT]we show the 
corresponding plots for the lines induced by the decay ^dm — > fi + fi~v, in the cases of an 
intermediate scalar and an intermediate vector particle, respectively. In these scenarios 
the expected line signal is somewhat weaker. 

The case mjv ~ m^, DM . Let us now turn to the case where the mass of the 
neutral fermion produced in the tree- level decay i/jdm — > £ + £~N is comparable in size to 
the dark matter mass itself. This possibility can occur, for example, within the leptophilic 
model discussed in section [2. 2. 2\ where -^dm is the hidden gaugino of an unbroken £7(1)- 
symmetry, and A is a neutralino [9] . 

As was shown in section 12.2. 1[ the decay channel ^dm - ► 7 A is kinematically en- 
hanced compared to the three-body decay when ~ m^, DM , provided that V'dm and 
N have opposite CP parities (77 = —1). Thus, such scenarios can be tested particularly 
well via the loop-induced gamma-ray line signal. In order to infer the observational con- 
straints, it is convenient to consider the ratio m^ DM /r(?/> DM — > jN), which determines 
the magnitude of the observable flux. For the case of scalar-mediated decay, and purely 
chiral couplings, it is given by (see eq. (12.311) ) 

[ ; feM i? chir V 10 26 sx2.5TeV J 

{ 3 x 10 29 s TeV for m N -> 0, 77 = ±1 

7xl0 29 sTeV for m N -> m^ DM , rj = +1 ; (4.9) 

1.4 x 10 28 sTeV (^-) 2 for m N m^ DM , 77 = -1 
V +DM / 

where R chlT = 1 for three-body decays into a single lepton flavor, and i? chir = 3 for flavor- 
democratic three-body decays. In the case of an intermediate vector, the right-hand side 
is smaller by a factor nine, implying a nine times larger gamma-ray flux. From the last 
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line, it is apparent that in the case of opposite CP parities, the monochromatic gamma- 
ray flux is enhanced for large values of m^ DM , when keeping the photon energy E y fixed. 
Note that a similar enhancement of the decay channel ipBM — > Z°N, which may also be 
induced at the loop-level, could lead to complementary constraints from the antiproton 
flux produced by the fragmentation of the Z-boson, which we do not discuss here. 

In order to illustrate this result, we consider a number of benchmark scenarios for 
which m^ DM and are of comparable size, with parameters chosen as shown in table [TJ 
All the benchmark scenarios reproduce the PAMELA positron data, and all except sce- 
narios 1, 2 and 3 additionally reproduce the electron spectrum measured by Fermi. Note 
that the maximum lepton energy in the three-body decay ipr>M — > i + £~N coincides with 
the energy of the monochromatic photons, E max = E 7 . 

The gamma-ray line signal induced by the one-loop decay ^dm — > 7 A is shown in 
fig. [12] for the various benchmark scenarios. Clearly, the scenarios 1 and 4 are in conflict 
with the gamma-ray line searches performed by Fermi and HEGRA, respectively. Thus, 
despite the fact that the dark matter couples only to leptons at tree-level, the gamma-ray 
line signal induced by one-loop corrections has an intensity that is detectable by present 
gamma-ray telescopes. In other words, scenarios 1 and 4 can be ruled out as possible 
explanations of the high-energy positron excess, because the loop-induced radiative decay 
produces a gamma-ray line that should have been already detected. This shows that 
the higher-order corrections are indeed relevant and have to be taken into account. In 
contrast, the other benchmarks are in agreement with present bounds on gamma-ray 
lines. For example, in scenario 3 the partial lifetime for the radiative decay is larger 
compared to scenario 1, and lies slightly above the current Fermi bounds. Scenario 6 
can be tested in the future by the CTA. Since there is no kinematic enhancement of the 
decay V'dm — > 7 A" in the case 77 = + 1 , the intensity of the gamma-ray line is comparably 
weak in scenarios 2 and 5. For example, scenario 5 differs from 6 just by the sign of rj, but 
is much more difficult to probe by the CTA. Finally, for benchmark point 7, we assume 
that the couplings of the leptons to -^dm and to N have opposite chirality, in which case 
the loop is strongly suppressed and there is no hope of detecting a gamma-ray line signal. 



5 Conclusions 

We have analyzed the radiative decay of dark matter particles in view of the leptonic 
cosmic-ray anomalies reported by PAMELA and Fermi LAT. Assuming an effective de- 
scription of leptophilic dark matter decay, we have pointed out that the lines induced 
at the quantum level may be observable and can be used to constrain models of de- 
caying dark matter. In the case of scalar dark matter, two-body decays into photons 
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are strongly helicity-suppressed and thus unobservable. In the case that the dark mat- 
ter particles carry spin 1/2, however, the radiative decay rate is typically suppressed 
compared to the tree-level decays by some two to three orders of magnitude. Interest- 
ingly, the corresponding partial lifetimes for decays into monochromatic photons can 
then be in the observable range, in particular for dark matter masses of a few hundred 
GeV, where stringent constraints from Fermi LAT apply. Thus, in some cases the loop- 
induced gamma-ray line yields constraints that can be competitive with the constraints 
on charged cosmic rays. At higher energies, constraints from Cherenkov telescopes ex- 
ist. At present, these constraints are only relevant for certain scenarios for which the 
radiative two-body decay is kinematically enhanced compared to the three-body decay 
channel. However, we have pointed out that the proposed CTA should be able to improve 
on the existing bounds significantly and probe a relevant part of the parameter space 
which is presently unconstrained. 
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A Decay widths for fermionic dark matter 

A.l The decay ^ D m £ + £~N 

The differential decay rate for this process is given by 

dT(ij DU -)■ e+e-N) = 1 \ \M t + M u \ 2 dtds . (A.l) 

(2tt) 3 64m^ DM 

Note that there is a relative minus sign between the t- and u-channel amplitudes due to 
the exchange of two anticommuting fermions that is not present by a naive application 
of the Feynman rules for the two diagrams. Neglecting the lepton mass, one obtains for 
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the squared amplitude 



\M t + M u \ 



L |2 



+ \ A %?) {\ A w\ 2 + l^«vl 2 ) 



(t - m 2 N ){ml DM - t) (u-m 2 N )(ml DM -u) 



(t - m\Y 



(u — m|) 2 



2r]{ Re 



(\L*\L Y , (\R*\R \ 2 
\ A tip A W) \ A li> A lN) 



m^ DM m N s 



(t-m|)(u-m|) 



< (t - m 2 N )(m 2 buu -t) + (u- m^)(m^ M - n) - s(t + n) 
(£-m 2 )(u-m|) 



(A.2) 



where 



(ft- Pi) 2 , £=(<?i- P2) 2 , u = (qi -P3) 2 = ml + m 2 N + 2m 2 - s - t . (A.3) 



Again, 77 = r]^ DM r] N = ±1 depending on the CP eigenvalues of -0dm and N. The integra- 
tion limits for the Mandelstam variables are given by 



< s < (m. 



<pDM 



m N f 



and 



where 



L > z 2 V TOM 



+ m 2 N - s =f ^ A ( m J DM ' m &' s )) 



A(a, 6, c) = a 2 + 6 2 + c 2 - 2a6 - 2ac - 26c 



(A.4) 

(A.5) 
(A.6) 



We can perform the kinematical integrations in the limit 3> t, u, in which case the 
Mandelstam variables in the denominator can be neglected. We then get 



r(^DM i + i~N) 



1 



TOM 



64(2vr) 3 6m| 
- ?7Re (A^A^yA^A^r) 



(\ A ei>\ 2 + l^>| 2 ) (\ a £n\ 2 + l^wl 2 ) 



+ 2r/Re 



F 2 (x) 



(A.7) 



where x = m^/mjl and -Fi(x), -F 2 (a;) are defined in eqs. (12.61) . (12. 7p . 

In the case of mediation by a vector, the matrix element for vanishing lepton mass 
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reads 



\M t + M U \ 2 = 4 (|A,V^| 2 + |A* Xf N \ 2 ) 



x 



(u - m 2 N )(m 2 , om -u) | (t-m^)(m^ DM -t) 



(t - m 2 ) 2 



+ 4 ( I A^,A^v| + | A^A 



R \ L |2 

lip A IN I 



+ 8rW Re 



s(t + it 



(u — m v ) 2 
1 



+ 



(t — m v ) 2 {u — m v ) 2 
m ^BM m NS 



- (t - m 2 v )(u - m 2 v ) 



+ 2Re [A^A^Aj^A^] 



s(t + u) 



(t-m v ){u- 



m\ 



(A.I 



In the limit m^ DM <C my we get for the decay rate 

1 



64(2vr) 3 6m^ 
+ 277 Re (A^A^AjAf^) J 



(\ A £^\ 2 + |A^,| 2 ) (|A^v| 2 + |A^v| 2 ) 



(A.9) 



A. 2 The decay ^ D m — >■ lN 

There are four scalar-mediated diagrams at the one-loop level contributing to the de- 
cay ipr>M — > ^N, which are shown in fig. [2J Due to gauge invariance, in the case of 
CP-conserving interactions, the matrix element corresponding to the sum of the four 
diagrams can be written in the form 

M = l -^u(h)(P R - VNm P L )a^k 2 ,e:u(p) 

= -^^u(h)(P R - VNvMhfuip) , (A.10) 

where a^ v = i[y^, 7 v ]/2 and % DM , r] N are the CP eigenvalues of V'dm and N, respectively. 
This is manifestly gauge invariant in the sense that it satisfies the Ward identity: the 
matrix element vanishes when replacing e* —> k 2fl since the photon is on-shell. 

The effective coupling g]v 7? /, DM for an intermediate scalar can be given in terms of 
loop integrals as follows, 

+ ( A w x w - ^DMWAfivA5,)[^ DM m^ DM (/ 2 — K) — r) N m N K}\ , (A. 11) 
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where the sum runs over all fermions / and all mediators £ that contribute in the loop. 
The loop integrals are written in terms of Feynman parameters as 



1 

A 



dx 
1 — x 



hgX 



K 



— / dx\ogX 



—— I dx I 1 + 



B 



A x(l — x) 



\okX 



(A.12) 
(A.13) 
(A.14) 



where 



A 
B 

X 



2 2 

m, h — m N 

TOM " 

mix + m^fl 



mjx + 771^(1 



X) 

■ x) 



X) 

■ x) 



In the limit m^ E 



mjx + — x) — m 2 N x(l — x) 
rax <C ?7i£, the loop integrals take on the simplified form [48] 



mi 



2m\ 



f 2 (m 2 e /ml), 



k= 1 -p 

2 

where the functions /, f 2 are defined as 

f( x ) = zr L- i + 
h{x) 



i 



1 — X 

1 

(l-x) 



I — X 

1 + x4 



ln(x) 
2z 



x 



\n(x) 



(A.15) 
(A.16) 

(A.17) 



(A.18) 
(A.19) 
(A.20) 

(A.21) 
(A.22) 



The expression for the effective coupling then assumes the form 



V 



16tt 2 



f{m\jm\ 



mi 



Ft, \R \ ^DM™TC 



rj N m N 



4m| 



f 2 (m 2 :/m 2 



ey 2 

g4 vr 2 m V'DM 



where in the last line we have taken m^ — > 0. 
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f 1 - £ ^ ( (AfcAfc - . A&Aft) } , (A.23) 

V TDM / , E '' 4 S 



In the case of an intermediate vector, the same integrals /, J 2 , K (with the replace- 
ment mj — > my in the constants B, X) appear, together with the additional integral 



J 



1 f dx 



A x 



logX 



which simplifies in the limit m^ DM ,m N <C m s to 

j_ 1 ln(x) 7 _ 1 



m v 1 — x 



mi 



-f{m\jm 2 v 



where in this case the kinematical functions are defined as 



fix) 

Furthermore, one finds 



1 — x 
1 

(l-x) 



1 + 



x 



ln(x) 



1 — x 

5x 2x(l - 2x) 



1 3 + 3(1 - x 



ln(x) 



J-K =K- J 



Amy 



The effective coupling in terms of loop integrals is given by 



V _ e lN^ M 

9n^ — 



8^ 



+ ri N m N (J - K)\ + 2m < (7/0 DM A^ v Ag, - r] N \f N X^)J 
For m^ DM <C my this expression then simplifies to 



v 

9N-yip — 



e ^m^ MSr f ( L R R L ,f v {m,/m 2 v ) 



mi 



(A.24) 



(A.25) 

(A.26) 
(A.27) 



(A.28) 



(A.29) 



3 (vNV^oM^W^hp ~ A^atA^,) 



VNrriN 



Amy 



fUrn 2 /m 2 y)} 



3er] 



m 



fen 



32tt 2 m 2 . 



r\m^ 



where in the last line we have taken — > 0. 

The decay rate in both cases is finally given by 



(A.30) 



r(^DM -> 7^0 



57T 



A' 



m 



(A.31) 



'/'DM 
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For the scalar in the case of one mediator coupled to leptons we get in the limit mi <Cmjv 
and m^ DM < m s , 

3 



r(v>DM jN) 



m 



♦dm 



87r(647r 2 ) 2 m| 



1 



N 



711 



1 



r]m N 



X 



whereas for the vector we have, in the limit Cmjv and ?n,/, DM <C ^y, 



1 9e 2 m5 pM 
Stt (8vr 2 ) 2 16m^ \ m 



N 



♦dm 



1 - 



r]m N 



y^ (A^yA^, ?7 A^A^,) 



(A.32) 



(A.33) 



B Decay widths for scalar dark matter 



In this appendix we present the expressions for the decay width of the radiative decay 
of scalar dark matter into two photons, 0dm — 1 77- 

For Xf, = Xf^ = Xty the decay rate reads [50j EE] 

2 



'DM 



-> 77) 



m 



TOM 



,2 \ 2 



4tt V 16tt 2 
where = ^| DM /(4m 2 ) and 

A / (r) = 2[r + (r-l)/(r)]/r 2 

/(r) = 



arcsm Wr 



In V , — Z7T 



r < 1 



t > 1 



i-v/T^V 

In the case of interest here, rg 3> 1. Then we can approximate 

>2 



A/(r) ~ - \ 2 - ^(ln(4r) - ztt) 



fB.r 



(B.2) 
(B.3) 



(B-4) 



t y 2 

In this limit, and taking only one lepton species into account, the decay rate is given by 



-> 77) 



|A^| 



,2 \ 2 



4m» 



167T ■'" WJM V16tt 2 / m 2 

<PDM 



7T 2 1,9/ 

2 + y-2 ln ( 4 ^) 



vr 2 ln 2 (4r,) 



(B.5) 
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